Astronomy & Astrophysics manuscript no. 


February 1, 2008 


(DOI: will be inserted by hand later) 





Cosmic Background dipole measurements with Planck-High 

Frequency Instrument 

M. Piat\ G. Lagache\ J.P. Bernard\ M. Giard^, and J.L. Puget^ 

^ Institut d'Astrophysique Spatiale (IAS), Universite Paris Sud, Bat. 121, 91405 Orsay Cedex, FRANCE 
^ Centre d'Etude Spatiale des Rayonnements (CESR), 9 avenue du Colonel Roche BP 4346, 31028 Toulouse 
Cedex4, FRANCE 

Received October 19, 2001; accepted ??? ??, 200? 

Abstract. This paper discusses the Cosmic Background (CB) dipoles observations in the framework of the Planck 
mission. Dipoles observations can be used in three ways: (i) It gives a measurement of the peculiar velocity of our 
Galaxy which is an important observation in large scale structures formation model, (ii) Measuring the dipole can 
give unprecedent information on the monopole (that can be in some cases hard to obtain due to large foreground 
contaminations), (iii) The dipole can be an ideal absolute calibrator, easily detectable in cosmological experiments. 
Following the last two objectives, the main goal of the work presented here is twofold. First, we study the accuracy 
of the Planck-HFI calibration using the Cosmic Microwave Background (CMB) dipole measured by COBE as well 
as the Earth orbital motion dipole. We show that we can reach for HFI, a relative calibration between rings of 
about 1% and an absolute calibration better than 0.4% for the CMB channels (in the end, the absolute calibration 
will be limited by the uncertainties on the CMB temperature). We also show that Planck will be able to measure 
the CMB dipole direction at better than 1.7 arcmin and improve on the amplitude. Second, we investigate the 
detection of the Cosmic Far-Infrared Background (FIRB) dipole. Measuring this dipole could give a new and 
independent determination of the FIRB for which a direct determination is quite difficult due to Galactic dust 
emission contamination. We show that such a detection would require a Galactic dust emission removal at better 
than 1%, which will be very hard to achieve. 



1. Introduction: the Cosmic Background and its 
dipole signal 

The Cosmic Background (CB) is the extragalactic part of 
the diffuse electromagnetic emission at all wavelengths. If 
. the universe obeys the cosmological principle and is ho- 
' mogeneous and isotropic, this background is expected to 
be nearly isotropic in the rest frame where the matter in 
a large volume around the observer has no bulk velocity 
other than the Hubble expansion. The measurement of the 
intensity and Spectral Energy Distribution (SED) of this 
isotropic background is a difficult observational challenge. 
It requires a separation of this extragalactic component 
from all other diffuse foregrounds (interplanetary or in- 
terstellar emission). Although the CB was detected quite 
early in the radio (radio-galaxies), Xrays and gamma rays 
as it dominates over the foregrounds, these components 
account for only 0.027% of the electromagnetic content of 
the present universe. The CB is dominated by the Cosmic 
Microwave Background (CMB) which accounts for 93% of 
electromagnetic content of the present universe and is a 
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truly diffuse component coming from the pregalactic era 
of the universe. The second component in energy content 
is the radiation from galaxies integrated over all redshifts 
in the ultraviolet-optical-near-infrared (stellar emission) 
and in the thermal and far infrared (re-emission of stellar 
radiation absorbed by interstellar dust). This component 
has been detected only recently, first in the far infrared 
in the COBE-FIRAS and COBE-DIRBE data (Puget et 
al. 1996; Fixsen et al. 1998; Hauser et al. 1998; Lagache 
et al. 2000) then in the optical with HST (Bernstein et 
al. (in preparation, quoted by Madau & Pozzetti 2000); 
Pozzetti et al. 1998) and near infrared (Gorjian at al. 
2000; Dwek & Arendt 1999; Cambresy et al. 2000; Wright 
2000). These measurements have been complemented by 
very high-energy gamma rays absorption by this CB (see 
for example Renault et al. 2001) but are still affected by 
rather large uncertainties especially in the optical and mid 
infrared. Figure |l| shows the spectrum of the CB from ra- 
dio to gamma rays (data from Gispert et al. 2000 for the 
part [30nm; 1mm] and from Halpern & Scott 1999 for the 
radio part). 

The CB should display a dipole anisotropy associated with 
the peculiar motion of the observer with respect to the lo- 
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cal cosmological standard of rest. This peculiar motion has 
been first detected for the CMB by Kogut et al. (1993) . 
It is composed of several terms: the peculiar motion of 
the Sun (sum of the Galaxy peculiar velocity and the ve- 
locity of the Sun with the rotation of our Galaxy), the 
Earth velocity in its orbit around the Sun and the spe- 
cific motion of the observer with respect to the Earth 
(orbital motion of the satellite for example). The main 
term is the first one with about 369 km/s in the direc- 
tion (/,6) ~ (264.3°, 48.0°) (Lineweaver et al. |1996| ). The 
second one, dominated by the Earth orbital motion, is a 
sinusoidal term with amplitude of about 37 km/s . It is 
thus about 10% of the first one but has a very specific 
character it is perfectly known in amplitude and direction 
and changes periodically with time during the year. The 
COBE-DMR experiment has measured these terms accu- 
rately. Fixsen et al. (1996) has shown that the CMB spec- 
trum absolutely measured with the COBE-FIRAS experi- 
ment with high accuracy could also be measured using its 
dipole anisotropy. The two spectra obtained were compat- 
ible within error bars. The peculiar motion of the galaxy 
has also been measured with respect to distant galaxies as 
an an isotropy in the Hubble constant (Lauer & Postman 
1994). Although the accuracy is not very good, the value 



obtained is compatible with the one measured by COBE 
(Scaramella et al. 1991, Kogut et al. 1993). These measure- 
ments confirm that the rest frame defined as the one in 
which the CMB is isotropic coincides with the rest frame 
defined by the galaxies on large scales. 
Motion with velocity (3 = v/c through an isotropic radi- 
ation field of intensity /(j/) yields an observed intensity 
given by: 



13 cose)] 
(l-/3cos6')3 



(1) 



where v is the frequency and 9 is the angle between (3 and 
the direction of observation, these two parameters beeing 
measured in the observer's frame. The spectral intensity 
of the dipole amplitude is therefore given by: 

/dip(j^) = /obs(t^, = 0°)- /obs(j^, e = 90°) 

~ 13 cosd (^I{v) - ^'^{^)^ for /3 < 1 (2) 

The spectral intensity of the CB dipole amplitude is shown 
in figure |l] and its relative level in figure ||. 
The dipole anisotropy of the CB has a two-fold interest 
for future CB projects: 

— At frequencies where the separation of the CB from 
foregrounds is difficult (from submillimeter to ultravio- 
let), an accurate measurement of the dipole anisotropy 
will give the SED of the CB independently of the in- 
terstellar and interplanetary foregrounds which do not 
show a dipole anisotropy. 

— A very accurate measurement of the CMB dipole 
anisotropy over a one year period should allow to use 
the dipole anisotropy associated with the orbital mo- 
tion of the Earth around the Sun as the best photomet- 
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Fig. 1. Spectrum of the CB (top) and its dipole amplitude 
(bottom) from radio to ultra-violet. 
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Fig. 2. Amplitude of the CB dipole relative to the CB 
intensity expressed in % from radio to ultra-violet. 

ric calibrator for extended sources in the microwave re- 
gion of the spectrum. It should therefore improve the 
determination of the peculiar motion of the solar sys- 
tem. This is possible due to the Planckian nature of 
the CMB and the accuracy of its temperature mea- 
sured with the FIRAS spectrometer. The CMB pro- 
vides today an extended photometric standard with an 
accuracy of less than 2 10^^ in thermodynamic tem- 
perature (Fixsen et al. 1996). 

This paper explores these two questions in the context 
of the Planck mission and specifically its High Frequency 
Instrument (Planck-HFI, Lamarre 199S). It is divided in 
2 parts. The first one deals with the calibration of Planck- 
HFI using the CMB dipole. In section 2, we present the 
calibration philosophy of Planck, while section 3 concerns 
more specifically the calibration accuracy. The second 
part, in section 4, is a study of the detectability of the 
Far InfraRed Background (FIRB) dipole. 
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2. Planck-HFI calibration philosophy 

Different sources will be used to monitor inflight the re- 
sponse or to get the photometric calibration of the Planck- 
HFI channels: 

— Extended sources: the CMB dipole, the Earth orbital 
motion dipole, the Galactic disc and eventually maps 
of Herschel-SPIRE. 

— Point sources: planets, infrared galaxies and asteroids. 

A general principle throughout the calibration for all chan- 
nels is that relative calibrations will be established as pre- 
cisely as possible between individual scans before estab- 
lishing an absolute calibration for all data. 
For channels dominated by the CMB (lOOGHz to 
217GHz), as will be seen in next section, the dipole sig- 
nal from the CMB will be used to measure the relative 
response variations on all useful timescales since it is an 
extended source available all over the mission. The ampli- 
tude and direction of the cosmic dipole has been measured 
by COBE. It will be remeasured independantly by Planck 
using the dipole from the Earth orbital motion as the pri- 
mary calibrator of the instrument absolute response. 
For higher frequency channels, mainly dominated by 
galactic dust emission, the calibration will be done by 
comparison with COBE-FIRAS map. The absolute cali- 
bration will also be done on the whole sky. The relative 
variations on shorter time scales will be measured inde- 
pendantly and are out of the scope of this paper. A short 
summary of the technique is given below. 
The response variation of a bolometer depends on a few 
fixed (or very slowly varying) parameters like the thermal 
conductance or the bias current for exemplc. It also de- 
pends on changing parameters, mainly the temperature 
of the heat sink of the bolometer and the incident power 
background. This last quantity is directly related to the 
flux from the sky which is known to first order before 
the mission or after a first iteration on the data, and to 
the temperature of the various stages of the optics (1.6 
K, 4K and the telescope). Since these temperatures are 
monitored during the mission, it will be possible to corre- 
late the response measured as a function of time with the 
temperature of the relevant elements and adjust very pre- 
cisely the parameter of the model describing the bolome- 
ter chain. This can be done very accurately using the slow 
variations of the different temperatures. This model with 
the temperature measurements (or any other monitored 
parameter found to affect the response) can be used to 
interpolate, if needed, the short time scale response vari- 
ations on high frequency channels. Relative variations of 
the response will thus be established on many time scales. 
A preliminary absolute calibration will be built after typ- 
ically a month period and will be improved until the final 
full data reduction for the detector is done. The bolome- 
ter model and the temperature measurements will be used 
to interpolate for time scales between one hour and one 
month. 

The CMB dipole is a strong known source (using COBE 



data initially or its properties remeasured by Planck at the 
end of the mission) which can be used to monitor the re- 
sponse variations over various time scales for all channels 
where its signal either dominates or is clearly detectable. 
It is thus important to evaluate the achievable accuracies 
for both the final absolute calibration and the measure- 
ments of the relative response. 

3. Using the CB dipole as a calibration signal for 
HFI 

The CMB shows anisotropics on all scales tracing the 
small inhomogeneities in the pre-galactic era that lead to 
the formation of structure in the universe. These fluctua- 
tions are observed on a surface of last scattering as seen 
from the Earth. They contain high spatial frequency terms 
which are the main objective of the Planck mission and 
low frequency components which arc not carrying much 
cosmological information as they are much more variable 
with the observer position (cosmic variance). The dipole 
term of this cosmological component is indistinguishable 
from the one associated with the peculiar motion of the 
solar system and is therefore a fundamental limit for the 
accuracy to measure this motion. The CMB cosmologi- 
calanisotropies are thus defined as with no dipole compo- 
nent (^ = 1 term of the usual decomposition in spherical 
harmonics of the sky) . The contribution of higher I to the 
dipole measurement is zero when it is done on the whole 
sky (by definition) but is not zero when the dipole is fitted 
on a fraction of the sky. In this section, we first present 
the Planck-HFI simulations we used, and explore the rel- 
ative calibration on rings as well as the global absolute 
calibration accuracy. 

3.1. Simulations 

The Planck satellite will orbit around the lagrangian L2 
point and will therefore follow the Earth orbital motion 
around the Sun. Its scanning strategy consists in rotat- 
ing the satellite at 1 rpm around its spin axis which will 
follow the antisolar direction. The beam axis, located 85° 
away from the spin axis, scans the sky along large circles. 
We define a "circle" as beeing one instantaneous Planck 
observation of 1 minute. The spin axis is relocated about 
every hour so that each circle is scanned about 60 times 
covering the whole sky on six months. These 60 circles 
are averaged into what we call one "ring". The spin axis 
follows a sinusoid trajectory along the ecliptic plane with 
6 oscillations per year and with full amplitude of ±10° so 
that the polar caps are not left unobserved. 
For simplicity and to spare computer time and place, only 
84 rings uniformly distributed over the sky in a one year 
mission were simulated. Moreover, the removal of slow 
drifts on circles has to be considered together with the 
process of removing systematic effects using redundancies 
and is therefore out of the scope of this paper. 
The simulated sky is the sum of 5 contributions: 
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1. CMB cosmological anisotropies, obtained from a stan- 
dard Cold Dark Matter model: fltot = 1, fib = 0.05, 
r^A = and Hf) = 50km/s/Mpc. The Ce power spec- 
trum is computed with CMBFAST (Seljak et al. 1996) 
and the map is generated in a healpix-type all-sky pix- 
elisation with a pixel size of 3.5 arcmin (Gorski et al. 

1998) . 

2. The solar system Peculiar Motion dipole (PM dipole) 
assuming an amplitude of 3.36 mK in the direction 
{l,b) = (264.31°, 48.05°) (Lineweaver et al. ^996|). 

3. The Earth Orbital Motion dipole (OM dipole) assum- 
ing an amplitude of about 336 /iK with its maximum 
at zero ecliptic latitude in the direction of the OM. 

4. Galactic dust, scaled from IRAS 100/im map using 
a 17. 5K blackbody modified by a ly^ emissivity law 
(Boulanger et al. 1996). 

5. One realisation of the noise with levels given in table 
|l|. These figures were estimated from the sensitivities 
of t he Pl anck-HFI channels as defined by Lamarre et 
al. ( |2001 ). We include a 1//^ contribution to the noise 
power spectrum with a knee frequency of lOmHz. In 
order to define the noise level for rings, the instanta- 
neous sensitivity on a circle has been divided by 
assuming that all circles are coadded into a single ring 
and that the noise is not correlated from one circle to 
the other. This is not strictly the case given the fre- 
quency aliasing produced by the co-addition (Janssen 
et al. |199(: , Delabrouille 1998b| ). However, it is a rather 
good approximation (see figure 3 from Giard et al. 

1999) . Slow drifts at frequencies lower than the spin 
frequency /spin = l/60Hz are assumed to be removed 
by a destriping algorithm (Delabrouille |1998c| ). 

3.2. Ring Analysis 

For each rings, a linear regression is done between the 
simulated signal and the expected dipole signal which is 
the sum of the PM dipole and the OM dipole as shown 
on figure The slope of this linear regression gives the 
absolute calibration factor with its uncertainty. The varia- 
tions of this factor gives the relative calibration accuracy. 
The standard deviation of the calibration factor is there- 
fore a good estimator of the relative calibration accuracy. 
The absolute calibration will be determined finally on the 
whole sky map (see section ^3] for a first attempt). A mask 
on region with IRAS 100/xm emission higher than a cut-off 
value can be applied in order to decrease the contribution 
from dust emission. 



3.2.1. Relative calibration accuracy 

In this analysis, we removed regions with IRAS lOO/im 
emission higher than lOMJy/sr which is quite optimal 
as will be shown in 3.2.2. Figure ^ shows the varia- 
tions of the calibration factor for a HFI lOOGHz chan- 
nel along the simulated mission assuming known OM and 
PM dipoles characteristics, with and without CMB cos- 



mological anisotropies. While the PM dipole varies along 
the mission, the OM dipole produces almost the same sig- 
nal on each rings, as shown on figure ^. The accuracy with 
which the dipole can be measured remains high even when 
the PM dipole is at its minimum, on rings around num- 
bers 40 and 80. When CMB cosmological anisotropies are 
not included, the relative variations of the calibration fac- 
tor are very low, mainly dominated by noise effects. CMB 
anisotropies dominate the uncertainties in measuring the 
amplitude of the dipole on individual ring. The accuracy 
improves with averages on larger number of rings. This 
part of the error goes exactly to zero when the whole sky is 
used. The accuracy of dipole fitting on HFI rings is there- 
fore highly limited by the CMB cosmological anisotropies. 
This fit is equivalent to find the m=l mode of the ring 
power s pectrum Tm- The relation from Delabrouille et al. 
(1998a) links this spectrum with the Ce spectrum: 



E 



(3) 



^1 



where Wim are window functions related to the beam 
shape and to the ring radius. The dipole term given by 
Fi clearly depends on all multipole > 1 of the CMB 
cosmological anisotropies. This coupling between modes 
disappear on a full map, indicating that the final absolute 
calibration should be done on the full map of Planck data 



as illustrated in section 3.3. Nevertheless, relative calibra- 



tion between rings at lOOGHz is expected to be possible at 
a level of the order of a few percent as discussed in 3.2.3| . 
On 545GHz calibration data shown on figure H, the signal 
is dominated by galactic dust emission and the accuracy 
of the CMB dipole calibration is about the same consid- 
ering or not CMB cosmological anisotropies. The relative 
calibration between rings at 545GHz is possible at a level 
of about 15%. 

3.2.2. Shifted direction of the CMB dipole 

The PM velocity vector direction is known with an accu- 



racy of about 14 arcmin at la (Lineweaver et al. 1996). To 
test the effect of this uncertainty on rings relative calibra- 
tion, we fit the simulated data with a shifted PM dipole 
by 30 arcmin in all directions (about 2cr from COBE re- 
sults) . Figure ^ shows the results obtain close to the worst 
case. The accuracy on the relative calibration is not highly 
affected and the calibration factor standard deviation is 
increased by only 15%. These results show that a rela- 
tive calibration on rings is still possible in an iterative 
way: with a first step relative response correction and with 
enough observation time, the dipole direction will be im- 
proved and therefore the relative calibration will be en- 
hanced. 

3.2.3. Effect of dust contamination 

The effect of dust contamination on calibration can be 
tested by removing from the dipole fit the data points 
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Table 1. Expected sensitivities of Planck-HFI per detection chain. The last two lines give the expected sensitivities 
per channel for the full mission. 



Frequency (GHz) 


857 


545 


353 


217 


143 


100 


Resolution (arcmin) 


5 


5 


5 


5 


7.1 


9.2 


Number of detector 


4 


4 


4 


4 


4 


4 


NETcMB (/iK.Hz-"-^) 


182000 


3995 


553 


182 


123 


99 


Thermo, temperature sensitivity per ring (/iK) 


199000 


4370 


605 


200 


113 


80 


NEI (MJy.sr-\Hz-"-'') 


269 lO-" 


232 lO-" 


165 10"^ 


88 10-" 


4710"'' 


23 10"'' 


Intensity sensitivity per ring (MJy.sr"^) 


294 10"^ 


253 10"^ 


180 10~^ 


96 10"^ 


43 10"^ 


1910"^ 


Thermo, temperature sensitivity full mission (/iK) 


36500 


801 


111 


37 


17 


11 


Intensity sensitivity full mission (MJy.sr^^) 


5410"^ 


46 10"^ 


33 10^^ 


18 10'^ 


6.610"^ 


2.610"^ 
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Fig. 3. Example of the analysis realised on two rings. The black solid curve is the simulated signal. The OM and PM 
dipoles are represented by dotted and dashed curves respectively. The grey curve is the OM plus PM motion dipole 
fitted on the signal. 
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Fig. 4. Variations of the calibration factor for a HFI 
lOOGHz channel as a function of the simulated ring num- 
ber. We have assumed perfectly known OM and PM 
dipoles characteristics. The lower curve has been com- 
puted without CMB cosmological anisotropics and is 
shifted by 10% for clarity. Regions where IRAS lOO/im 
emission greater than lOMJy/sr were removed from the 
fit. 



Fig. 5. Variations of the calibration factor for a HFI 
545GHz channel as a function of the ring number. The 
legend is similar to figure ^. The lower curve is shifetd by 
200% for clarity. 



having an excess emission in IRAS 100/im map. For very 
low cut-off values, the fit is degraded due to the decrease 
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Fig. 6. Variations of the calibration factor for a HFI 
lOOGHz channel as a function of the ring number with 
an error in the PM dipole direction of 30 arcmin. The 
dotted curve assumed a perfect known PM dipole direc- 
tion. Regions where IRAS 100/im emission greater than 
lOMJy/sr were removed from the fit. 

in the number of data points. On the other hand, with a 
high cut-off value, all the emission of dust in our Galaxy 
is taken into account in the fit and it degrades the final 
accuracy. We therefore expect an optimum value to exist 
in between. In order to study this effect, the standard de- 
viation of the calibration factor on all rings is plotted in 
figure ^ as a function of the IRAS lOO/im cut-off. 
On CMB channels (lOOGHz, 143GHz and 217GHz), in- 
creasing the IRAS 100 fini cut-off value do not decrease 
the relative calibration accuracy: as shown in previous sec- 
tion, this accuracy is limited only by CMB cosmological 
anisotropics. A relative calibration accuracy of 1.2% at Icr 
is possible on these channels. On higher frequency chan- 
nels, the emission of galactic dust is more dominant and 
there is, as expected, an optimum cut-off value at a few 
MJy/sr. 

3.3. Absolute calibration accuracy on all simulated 
data 

In order to test the accuracy of absolute calibration on 
a global map, we have fitted simultaneously all simulated 
rings with a dipole term from OM and PM. This is an 
approximation of a full fit on the whole sky. Our procedure 
nevertheless gives a first approach to the global calibration 
which will be expensive in computer resources. 
We can write the signal as: 

Signal(?,&) = X [PM(^„,i„,,„)a,6)-|-0M(l,b) (4) 
+Ga\{l,b) + CMBA{l,b)] (5) 

where F is the calibration factor, Aq is the PM dipole 
amplitude {Ao—3.36 mK), {lo,bo) is the PM dipole direc- 
tion (^0 = 264.31°, 5o = 48.05°), Gal is the Galactic sig- 
nal and CMBA the CMB anisotropy contribution. This 
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Fig. 7. Variations of the standard deviation of the cali- 
bration factor (expressed in percentage of the calibration 
factor) for all HFI channels, except the 857GIIz, as a func- 
tion of the IRAS 100/im cut-off value. 

analysis has been done only on CMB channels (lOOGHz, 
143GHz and 217GHz). In order to remove dust contam- 
ination effect, we choose an IRAS lOOfim cut-off value 
of 15MJy/sr. Morevover, we remove CMB cosmological 
anisotropics from the fit since, in our case, they intro- 
duce a £ =1 term which will not be present when using 
the whole sky. We assume in the procedure that we per- 
fectly know only the OM dipole so that the calibration is 
due to this term. We first search for the best F, Aq, Iq 
and bo and their accuracy by fitting simultaneously the 
84 rings signal of one CMB channel. This analysis shows 
that we can separate our problem into two much faster 
fitting procedures: one to get the best calibration factor 
and PM dipole amplitude, and the other to find the best 
PM dipole direction and amplitude. This two problems 
are mostly decorelated since a change in the PM dipole 
direction cannot be compensated by a change in its am- 
plitude (or in the calibration factor). These two analysis 
lead to contour plot shown in figure || and ^. We obtain 
the following results for one lOOGHz detection chain: 

- We can recover the PM dipole direction at better than 
1.7 arcmin (Fig. |^). 

- We find Ao = (3.374 ± 0.007) mK at 95% CL and 

- F is found to be equal to 0.996 ±0.002 at 95% CL (Fig. 
§). 

This shows on the one hand that Planck-HFI will be able 
to make a new and accurate determination of the PM 
dipole direction. On the other hand, by combining the 
error on F and Aq and considering the OM dipole as an 
absolute calibrator, we see that we can make an HFI low 
frequency absolute calibration at better than 0.4%. This 
first attempt is moreover quite pessimistic since it takes 
into account only part of Planck-HFI one year data. 
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Fig. 8. Contour plot in the l/h plan of the obtained 
with the fit of the PM dipole direction, expressed in per- 
centage of confidence level, for a lOOGHz channel. We as- 
sumed known PM dipole amplitude and instrument re- 
sponse for this fit, and no CMB cosmological anisotropics 
(see text). The asterisk point represents the best fit while 
the diamond is the input dipole direction. Regions with 
IRAS 100/im emission higher than 15 MJy/sr were re- 
moved from the fit. 
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Fig. 9. Contour plot in the F/Aq plan of the obtained 
with the fit of the PM dipole amplitude and instrument 
response, expressed in percentage of confidence level, for a 
lOOGHz channel. We assumed a known dipole direction for 
this fit, and no CMB cosmological anisotropics (see text). 
The asterisk point represents the best fit while the dia- 
mond is the input simulation. Regions with IRAS lOO/im 
emission higher than 15 MJy/sr were removed from the 
fit. 



4. Detectability of the dipole effect in the Far 
InfraRed Background (FIRB) 

The FIRB is due to the integrated emission of distant 
redshifted galaxies. Its SED between 100/xm and 1mm is 



well represented by the following expression (Lagache et 
al. 1999): 



liv) = 8.8 10" 



(6) 



where vq = lOOcm-i, Tq = 13.6^ and P^iT) is the Planck 
function. This expression has been obtained by averag- 
ing over multiple clean region, almost homegenously dis- 
tributed on the sky, where dust contamination is expected 
to be negligible. We therefore assumed that it represents 
the SED of the FIRB at rest. The FIRB dipole spec- 
trum can be deduced from its monopole spectrum by using 
equation || and ||: 



/dip FIRB (^) - ( - 1.4 



\ exp X — I 



(7) 



where x = /ii'/(fcTo). 



4.1. Signal to Noise ratio 

In order to evaluate the expected sensitivity of HFI on the 
FIRB dipole signal, we first consider the case of instrumen- 
tal noise having a pure white spectrum. As for the CMB 
dipole, the FIRB dipole will be best detected on a full sky 
map where no aliasing from CMB anisotropics will occur. 
An estimation of the Signal-to-Noise Ratio (SNR) can be 
done on the £ — 1 multipole of the angular power spec- 
trum decomposition. The expected Ci for a dipole signal 
having an amplitude A is given by: 



(8) 



Knox (1997) and Tegmark (1997) have shown that the ef- 
fect of uniform instrumental noise can be accurately mod- 
eled as an additional random field on the sky, with an 
angular power spectrum given by: 



Ci 



(9) 



where ilj, is the beam solid angle and <j is the r.m.s. 
noise per pixel. The SNR on Ci is therefore equal to 
SNR = Ci/Ci noise and the results are summarised in 
table H It shows that the FIRB dipole can be detected at 
857GIIz and 545GIIz but not at 353GIIz. Even an averag- 
ing over 10° could allow to detect the FIRB dipole effects 
with a significative SNR on these two high frequency chan- 
nels if the detection noise is white. 

Nevertheless, the HFI will exhibit 1/f noise mainly due 
to the readout electronics (Gaertner et al. 1997, Piat 
et al. 1997), temperature fluctuations of the cryogenic 
system (Piat et al. 1999) and from far side lobe sig- 
nal (Delabrouille 1998b). The required knee frequency is 
lOmHz for HFI in order to have enough stability on a 
timescale of 1 minute corresponding to the spin period 
of the satellite (Piat et al. 1999). Destriping algorithm 
can remove an important fraction of low frequency drifts. 



Delabrouille (1998c) has shown that removing fluctuations 
on timescales larger than the spinning period of Planck 
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Table 2. FIRB dipole level in the Planck-HFI channels and SNR for three pixelisations assuming a white detector 
noise only. We also assumed a 1.17 year mission, square pixels with size given by FWHM and bolometer noise given 
on table 



Frequency (GHz) 


857 


545 


353 


FIRB dipole intensity amplitude (MJy/sr) 


1.6 lO-" 


3.410"* 


4.910"" 


SNR per pixel 


0.04 


0.01 


0.00 


SNR per pixel of 1° 


0.50 


0.12 


0.2 


SNR per pixel of 10° 


5.0 


1.2 


0.2 


CiPiRB {{MJy/srY.rad') 


3.510"'' 


1.710"' 


3.410"^ 


Ci noise {{MJy/sry.rad-') 


6.110-'-' 


4.6 10"'-' 


2.310"^ 


SNR full map on Ci 


565 


36 


1.5 



is possible thanks to redundancies. A conservative analy- 
sis could be made by assuming that all slow components 
of the noise, except frequencies lower than the spin fre- 
quency, are not filtered. Equation |^ allows to express the 
ring anisotropy power spectrum with the angular 
power spectrum on the sky. The window functions Wtm 



are given by (Delabrouille et al. 1998a): 



(10) 



where 6 is the ring angular radius, is the beam re- 
sponse function (assuming a symmetric beam) and Vtm 
are given by the following expression where Pim are the 
associated Legendre polynomials: 



2£+l\ (£-m)l 



Pirn (cos ( 



for m > 0, 



47r J + \ , 

= form<0. 

Equation || can be used to deduce the angular power spec- 
trum of the instrumental noise from the Noise Equivalent 
Power (NEP) spectrum. 

A mode m on a ring has an equivalent frequency band- 
width of \/(2Tspin) where Tspm is the spin period. 
Furthermore, the NEP has to be filtered by the beam in 
order to obtain the noise spectrum projected on the sky. 
The NEP is therefore related to the ring power spectrum 
of the noise by: 

NEP\f)B\f) = 2T,pi„r„ (11) 

where / = m/Tspin is the frequency. B{f) is the beam 
filter given by the following expression in the case of 
a gaussian beam shape of Full Width Half Maximum 
FWHM = abVSh[2: 



B{.f) = cxp 



1 / ^CbTspin 

^2 V sine 



(12) 



Equation H can therefore be expressed in terms of NEP: 



NEP^f ^m/T,pi, 



2T. 



spin 



(13) 



t=\m\ 

where the sum is done on all accessible > 



The value of 



is given by the highest frequency /max transmitted by the 



1 000 r 



100 : 



1000 



Fig. 10. Effect of 1// noise on the noise angular power 
spectrum noise normalised at high for a knee fre- 
quency of lOmHz. The solid curve correspond to a power 
a — 1 (see equation |l^ while the dotted, dashed and 
dot-dashed lines are for a = 2, 3 and 4 respectively. 

readout electronics, about lOOHz for HFI that leads to 
^max = T'spin/max/ siu 9 = 6000 for a ring angular radius 
e = 90°. 

Equation O can be viewed as a matrix multiplication: 



NEP^ = M X Cr. 



(14) 



where NEP and Cnoisc are vectors containing NEP{f = 
m/Tspm) and noise respectively. M is a triangular matrix 
since Tm depends only on Ci with £ > |m|, and the solu- 
tion is therefore quite easy to obtain. With this method, 
we get noise assumiug a noise spectrum of Planck-HFI 
given by: 



NEP^if) = NEP( 



1 



An 



/ 



(15) 



where NEPq is the NEP value in the white part of the 
noise, /knee is the knee frequency and a a constant. The 
effect of 1// noise on the angular power spectrum is shown 
on figure ^ for a knee frequency of lOmHz in 4 different 
cases of noise behaviour. The SNR on Ci will be degraded 
by a maximum factor of about 270 which nevertheless al- 
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Fig. 11. Dust dipole directions superposed to the CMB 
dipole. The crosses are points obtain with COBE-DIRBE 
data at lOO/im, 140/im and 240/im with increasing size 
of the symbol, assuming a cut in galactic latitude of 
±15°. These points have been extrapolated to HFI 350/im 
(857GHz) and 550/im (545GHz) channels represented by 
diamonds with increasing size. 



lows a detection at more than 2a at 857GHz if foregrounds 
would be negligible.. 

4.2. Estimation of ttie total dipole spectrum 

In order to evaluate the effect of calibration and compo- 
nent separation, the relative level of FIRE dipole with re- 
spect to other components has to be evaluated. The total 
dipole signal on Planck-HFI channels is produced by the 
CMB, the FIRE and also from the inhomeneously dust 
distribution in the Galaxy. 

The repartition of dust in the Galaxy produces a dipole 
signal, the so-called "dust dipole", that is not due to 
Doppler effect. This signal has been estimated on COEE- 
DIREE sky map at 100/im, 140/im and 240/im. In or- 
der to remove strong emission from the galactic plane, a 
cut in galactic latitude or a cut-off in maximum emission 
level has been applied. Both methods leads to about the 
same results for intermediate cut, within less than 20° on 
the fitted dipole direction and about 10% accuracy on its 
amplitude. The dipole that dust mimic is push toward 
the galactic south pole as the cut in the latitude is in- 
creased, or the cut-off in maximum emission is decreased. 
A reasonnable cut in galactic latitude of ±15° has there- 
fore been applied. These data have been used in order to 
interpolate the dust dipole on Planck-HFI 545GHz and 
857GHz as shown on figures ^ and |l2|. The dipole direc- 
tion has been extrapolated linearly with galactic coordi- 
nate, assuming a ±20° uncertainty. It seems to go in the 
direction of the south galactic pole, meaning that at higher 
wavelength only the north/south disymmetry will be seen. 
The dust dipole spectrum follows a (20 ± 1)K blackbody 
modified by a emissivity and with a H column density 




0.10 ^ 



0.01 



000 



Wavelength 'y/^rnj 



Fig. 12. Spectrum of the dust dipole. The diamonds are 
the points obtain with COEE-DIREE assuming a cut in 
galactic latitude of ±15°. The solid line is the best fit with 
a emissivity leading to a dust temperature of (20±1)K 
and a Hi column density of (5.2 ± 1.5) lO^^cm"^. Dashed 
lines gives the upper and lower limits. 




400 600 800 1000 

Frequency (GHz) 

Fig. 13. SED of dipoles from the CME (dotted line), dust 
(dipole projected on the CMB dipole direction, dashed 
lines giving lower and upper limits) and the FIRE (dot- 
dashed line). The solid region gives the spectrum of the 
total dipole with the uncertainty on the dust dipole, the 
black one beeing for positive dipole while the grey one is 
for negative dipole. 

of (5.2 ± 1.5) lO^^cm-^. We also assumed that 10% of this 
signal can be removed by component separation. The re- 
sulted dust dipole spectrum with its uncertainty, projected 
on the CMB dipole axis, is shown on figure |l^ together 
with the CMB, FIRE and total dipole. The FIRE dipole 
signal is only about 10% of the total dipole dominated by 
the dust emission one for frequencies higher than 700GHz, 
which has two consequences: 

1. To recover the FIRE dipole signal, the dust dipole has 
to be removed by componenent separation to a 1% 
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level accuracy (10 tinies better) which will be difficult 
to achieve, even on large scales. 
2. The calibration on the 545GHz and 857GHz chan- 
nels will be obtained by comparison of the galactic 
plane emission measured with COBE-FIRAS. The ac- 
curacy on these channels will probably be limited by 
the COBE-FIRAS calibration accuracy which is about 
3%, a factor 3 higher than needed to detect the FIRE 
dipole effect. 

The detection by Planck-HFI of the FIRB dipole is there- 
fore unfortunately very difficult, due to component sep- 
aration and calibration accuracies of the high frequency 
channels. 

5. Conclusion 

We have shown that we can use the PM and OM 
dipoles observation for relative and absolute calibration 
for Planck-HFI. On CMB channels, we reach a relative 
calibration accuracy of about 1.2% on rings and an abso- 
lute calibration accuracy better than 0.4% when the whole 
sky is used. Moreover, the CMB dipole direction will be 
improved by Planck. At the end, the Planck-HFI absolute 
calibration will be limited by the uncertainty on the CMB 
temperature. The next step consist of a detail study of 
systematics effects on calibration. 

The detectability of the FIRB dipole has also been stud- 
ied and we have shown that it will be a difficult challenge 
due to dust contamination. 
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